G3: Genes|Genomes|Genetics

Overexpression of -synuclein (-syn), an inherently disordered protein, triggers chronic activation of the mitochondrial unfolded protein response (UPRmt) pathway in Caenorhabditis elegans with enhanced dopaminergic (DAergic) neurodegeneration. Introduction of a loss-of-function(lf) mutation in atfs-1, the main transcriptional regulator of the UPRmt, into -syn nematodes results in significant neuroprotection from -syn-induced DA neuron loss, indicating that compensatory mechanisms provide neuroprotection. We performed a F3 forward genetic screen in C. elegans atfs-1(lf) mutants to identify molecular components associated with the modulation of neurodegeneration via UPRmt signaling in -syn-expressing DA neurons. Homozygous mutant animals were examined for enhanced neurodegeneration; multiple independent alleles were uncovered. Among these, we identified new nonsense alleles encoding the histone lysine demethylases (H3K27me3), jmjd-1.2 and jmjd-3.1. Another line carried a nonsense allele of twk-14. This gene encodes a conserved protein termed KCNK12 in mammals that facilitates passive background K+ leak currents to set and stabilize resting membrane potential. To further examine the association of these gene products in DA neurodegeneration, mutants and/or RNA interference were employed. DA neurodegeneration was observed in the -syn + atfs-1(lf) background when jmjd-1.2, jmjd-3.1, or twk-14 were individually depleted. These results provide evidence that jmjd-1.2 and jmjd-3.1, which encode previously characterized H3K27me3 demethylases, and the uncharacterized twk-14 gene product, orthologous to human KCNK12, naturally confer protection from -syn neurotoxicity.

Genetic-background studies require defined perturbations that can be crossed reproducibly into many recipient backgrounds. We generated a Drosophila dilp2GS-rpr donor line for adult-inducible ablation of insulin-producing cells (IPCs), which secrete insulin-like peptides and provide a tractable model of insulin-deficient metabolic physiology. This line carries dilp2-GeneSwitch-GAL4 and UAS-reaper in cis on the same second chromosome homolog over a balancer. PCR genotyping and sequencing confirmed both transgenic elements in the candidate recombinant line. RU486 induction reduced dilp2 mRNA expression, supporting partial IPC ablation. Treatment-duration testing identified 8 days of RU486 as sufficient to increase whole-body glucose in the dilp2GS-rpr line but not in the background-matched control; food intake did not differ between RU486- and vehicle-treated flies. Across metabolic assays, whole-body glucose showed the clearest RU486- and line-dependent phenotype. This validated dilp2GS-rpr line enables testing how recipient genetic backgrounds modify inducible IPC/DILP metabolic phenotypes and provides a framework for similar linked donor-line resources.

Pyruvate dehydrogenase complex (PDHc) deficiency (PDCD) is a primary mitochondrial disorder characterized by neurodevelopmental disability, altered intermediary metabolism and early mortality. Dichloroacetate (DCA), a pyruvate analogue, is a well-described PDHc activator that remains under clinical investigation for treatment of PDCD. Here, we studied the in vivo efficacy of a 5-point log concentration range of DCA on animal health and metabolism in C. elegans with feeding RNA interference (RNAi) expression knockdown of either PDHA-1 or DLD-1 homologues at graded degrees to model variable disease severity. These worm models recapitulate phenotypic features of PDCD observed in human patients, including reduced survival, delayed growth, locomotor impairment, and elevated lactate and/or pyruvate tissue levels. DCA treatment appeared well-tolerated, with no gross morphologic toxicity seen at doses up to 25 mM. Significantly improved health, survival, tissue lactate levels, and mitochondrial physiology were observed at 25 mM in pdha-1(RNAi) knockdown animals. DCA treatment in dld-1(RNAi) C. elegans models (undiluted, 1:20 dilution, and 1:100 dilution) showed significant therapeutic benefits on survival, neuromuscular function and metabolic phenotypes primarily in the moderate (1:20) and/or mild (1:100) dld-1(RNAi) deficiency strains, but not in full-dose dld-1(RNAi). Importantly, linear growth, neuromuscular activity, and mitochondrial physiology were significantly improved with DCA treatment even in the most severe dld-1(RNAi) undiluted model. Overall, preclinical modeling provides objective evidence of DCA therapeutic efficacy in C. elegans expression knockdown strains for two well-conserved homologues of PDHA1 and DLD that represent distinct genetic etiologies of PDHc deficiency, with demonstrated beneficial effects on survival, healthspan, tissue lactate, and mitochondrial physiology. These data further confirm that DCA's therapeutic effect correlates with PDHc disease phenotype severity in dld-1(RNAi) animals.

Characterizing the genetic basis of root system architecture and its role in early plant development is essential for developing maize varieties with improved nutrient uptake, enhanced early vigour, and higher yield potential in temperate regions. Landraces represent an invaluable source of allelic diversity that can be leveraged to enrich the genetic basis of modern breeding material. In this study, we used a high throughput phenotyping platform to characterize genetic variation for seedling root traits under chilling conditions relevant for early plant establishment in a large doubled haploid (DH) library derived from two European maize landraces. We dissected the quantitative genetic architecture of twelve seedling root traits using a haplotype-based genome-wide association study, identifying large-effect haplotypes specific to the individual landraces as well as numerous small-effect haplotypes present in both landraces. We validated the effects of four QTL in a biparental population, demonstrating their stability across genetic backgrounds. We found highly significant correlations between haplotype effects on seedling root traits evaluated in the phenotyping platform and early plant height evaluated in multi environment field trials, demonstrating the relevance of seedling root architecture for early plant establishment. In particular, haplotypes associated with seminal and lateral root length were the major determinants of early plant height under field conditions. Several of the haplotypes increasing seedling root length were absent from a broad panel of flint breeding lines, highlighting their potential as targets for introgression to improve early plant establishment under temperate growing conditions. Key messageSeedling root QTL discovered in a high throughput phenotyping platform under chilling conditions influence early plant development in the field.

The kynurenine pathway metabolizes tryptophan into 3-hydroxykynurenine (3-HK), a precursor for ommochrome eye pigments synthesized via the cardinal (cd) gene in mosquitoes. While cd disruption was presumed neutral, we observed fitness costs in Anopheles stephensi knock-in but not knock-out cd mutants. Here we investigated this anomaly further by assessing survival, fecundity, and midgut integrity across multiple cd mutant lines. Heterozygous knock-in lines, expressing a fluorescent marker and guide RNA for CRISPR/Cas9, exhibited reduced survival post-blood feeding, larva-to-adult survival deficits, and midgut barrier dysfunction, whereas knock-outs showed no such costs. Oral supplementation with xanthurenic acid partially rescued knock-in mortality, implicating oxidative stress linked to 3-HK metabolism. Expression analyses suggest transgene insertion effects, rather than cd disruption, underlie these fitness costs. These findings highlight the importance of evaluating insertional effects in gene drive target selection and support cd as a viable target for genetic control strategies in An. stephensi.

Hearing in Lepidoptera mediates a range of ecologically important behaviours, including mate communication, predator avoidance, and acoustic signalling. In moths, the evolution of predator-prey interactions with bats has further shaped hearing through a sensory arms race, with repeated co-option of auditory organs to detect and evade echolocating predators. Despite significant prior characterization of the neurophysiology and behaviour of hearing in moths, the genetic basis of hearing is poorly understood in most insects. In this study, we identify a core set of putative auditory genes in Lepidoptera using a combination of homology-based searches from Drosophila and evolutionary rate analyses. We find 56 genes present across all species and investigate whether gene copy number varies among non-hearing and hearing lineages and among 3 different ear types. We discovered seven genes associated with ear type and one with ear presence, but did not find significant losses in gene copy number in non-hearing species. We identified three genes (btv, Dnai2, and nompB) with strong evidence of selection in hearing clades and five genes with weaker evidence of selection. We discuss the potential roles of btv, nompB, and Dnai2 in ciliary transport and the aging of hair cells, as well as the possibility of actively amplified hearing. Our study serves as a primer and resource for further gene mining and functional testing of auditory genes in moths and other insects.

The P-element, a DNA transposon, has independently invaded two Drosophila species, accompanied by rapid evolution of suppression. In the germline, suppression is mediated primarily by maternally expressed piRNAs, a class of regulatory small RNAs associated with PIWI proteins. The offspring of females that lack P-element-specific piRNAs and males that contain P-elements suffer a syndrome of deleterious phenotypes, including sterility, genome rearrangements, gonadal atrophy, and mutations, while the offspring of the reciprocal cross are normal. These effects, collectively termed hybrid dysgenesis, have been investigated primarily in female D. melanogaster. Here, we study hybrid dysgenesis in male D. simulans. Using an attached-X chromosome stock, we generated genetically identical F1 males that differed only in maternal suppression of the P-element. Using targeted sequencing of P-element breakpoints, we show that P-element transposition is elevated in dysgenic males and confirm a preference for insertion near origins of replication. Using transcriptomics, we show that dysgenic males have elevated P-element expression and reduced splicing suppression, with patterns of gene expression suggesting the loss of mature sperm cells. Fertility assays show higher rates of male sterility but otherwise modest effects on fertility. In conjunction with the transcriptomic data, small RNA sequencing confirms that the piRNA pathway functions in testes. Our results suggest that the P-element may spread more readily through males than females, as transposition rates are similar while fertility defects are less severe in males.

The relationship between genotypic and phenotypic variation is determined by the complex interaction of genetic and environmental factors. While statistical methods capable of detecting such interactions exist, an axiomatic mathematical framework that seamlessly describes the combined effects of genetic modifications and environmental exposures on a common scale is lacking. In this report, buffering concepts are used to construct a measurement system that enables the geometric representation of both gene-by-gene and gene-by-environment interactions on the extended complex plane (i.e., as projections on the Riemann sphere). In this manner, any such interaction, or combination thereof, can be precisely defined and quantified as the deviation from the neutral value calculated through the applicable complex transformation. When thus conceptualized, the framework's parameterization defines the "state space" of a given measurable phenotype along both the real and imaginary dimensions, thus establishing an unambiguous and broadly applicable method for determining the phenotypic value expected upon combinatorial changes in genetic and/or environmental variables. Remarkably, by applying these methods, it is possible to quantify the effects of any gene-by-environment interaction using the equation, AGxE=Im([z]obs*zexp)/2, where zobs and zexp are complex numbers representing the observed and expected phenotypes of a given genotype expressed in terms of the buffering parameters, and b.

Septoria nodorum blotch (SNB) and tan spot (TS) wheat diseases are caused by necrotrophic fungal pathogens Parastagonospora nodorum (Pn) and Pyrenophora tritici-repentis (Ptr), respectively. Although recognised as premier model pathosystems for our understanding of necrotrophic effectors, no resistance mechanism has been reported in both diseases. Here, two SNB and TS resistance wheat lines ( 56:ZWB11 and 105:ZIF14) derived from the Australian national germplasm evaluation programme (CAIGE) were used to develop a double haploid mapping population. Two Pn and Ptr isolates of different pathotypes, their respective culture filtrates and effector SnTox267 were evaluated on the population. Genetic analysis of Ptr conidial inoculation of race 1 and race 2 identified a major resistance quantitative trait locus (QTL) (QTs.cur-1B) on chromosome 1B, while resistance to SNB was explained by several minor QTL. SnTox267 sensitivity was mapped to six locations (2A2, 2A3, 2B1, 2D3, 5B and 7B1) with only one QTL co-localized to known corresponding gene Snn7. Sensitivity loci 5B and 7B1 also conferred SNB resistance at seedling and adult stages. Two QTL on chromosome 2D1 and 7B2 were common in both SNB and TS, associated with disease at seedling stage and culture filtrate bioactivity, respectively. Resistance responses of 56:ZWB11 and 105:ZIF14 were confirmed cytologically, however, distinct responses were observed on wounded leaves. The defence responses were more effective against Ptr, while resistance to Pn infection was likely a combination of lack of susceptibility and effective physical barriers. Overall results demonstrated the distinction between the underlying resistance mechanisms to TS and SNB.

Survival analysis is a workhorse assay in Drosophila research to evaluate somatic fitness. It is indispensable in the study of ageing and insightful in immunity, metabolism, radiobiology, toxicology, ecology, and others. While conceptually simple, lifespan measurement is labour-intensive because it requires the continuous manual maintenance of large experimental cohorts. Here, we describe Drosben, an approach that combines a 3D-printed device to transfer flies from several vials simultaneously, a paper system for quick data recording and accompanying software that automatically digitalises life tables for analysis. We show that using Drosben reduces the time investment to perform lifespan assays by ~85%, with improved speed regardless of experience handling Drosophila vials. Using Drosben, we address the effects on longevity of chronic feeding of indole-acetic acid (IAA), naphthalene-acetic acid (NAA) and trimethoprim (TMP) -- compounds used to control heterologous targeted protein degradation systems. We find that IAA and NAA have noticeable deleterious effects while TMP has a small protective effect specifically in females. We further show that strong static magnetic fields do not affect Drosophila lifespan. Our work suggests that Drosben can cheaply accelerate research where lifespan is used as a life history trait.

Aedes aegypti is a major vector of arboviruses, and belongs to subfamily Culicinae, a diverse group of mosquitoes with homomorphic sex-determining chromosomes. Males are the heterogametic sex with a dominant male-determining locus (M locus). The M locus and its counterpart m locus are embedded in a region of suppressed recombination, with a large portion of this recombination desert showing significant molecular differentiation despite homomorphy. We developed a mathematical framework to examine M-linked genome editors that specifically target the m-chromosome during spermatogenesis, mimicking the naturally occurring sex ratio distorters (SRDs) in Culicinae that produce male-biased meiotic drives. Unlike previous models for species with heteromorphic sex chromosomes (e.g., X and Y), we incorporate features stemming from the homomorphic nature of the Ae. aegypti sex chromosomes such as varied linkage to the M locus, making the degree of super-Mendelian inheritance readily tunable. We evaluated in silico SRDs with a range of M-linkage and editing efficiencies and established the theoretical foundation for developing highly efficient SRDs that outperform several methods of population suppression. These SRDs can be tuned to mitigate impact on a neighboring population. The framework developed here is suitable for exploring SRD-mediated genetic biocontrol of pests with homomorphic sex chromosomes.

Horseweed (Erigeron canadensis) is a widely distributed annual weed that can cause significant yield losses if not properly controlled. Its phenotypic plasticity allows it to rapidly acclimate to new environmental conditions, such as drought and herbicides, such as glyphosate, with the potential for cross stress acclimatization. The objectives of this research were to uncover the physiological and genetic effects at the intersection of drought stress and glyphosate resistance. To this end, we performed greenhouse dose response experiments, RNAseq, 14C glyphosate absorption and translocation, and cuticular lipid profiling via GC/MS. Greenhouse dose-response experiments revealed that, after drought stress, there was a 2.5-3.7 fold reduction in glyphosate sensitivity via a significant reduction in glyphosate absorption, regardless if the starting population was resistant or susceptible to the field use rate already. Cuticular waxes were collected from each population with and without drought stress and were analyzed via GC/MS. When comparing total wax loads of plants grown under WW and DS conditions, we found that drought stress significantly increased total wax loads for all three populations. Additionally drought stress substantial increases the proportion of triterpenoids in the cuticle. By RNAseq, we found serval triterpenoid biosynthesis genes upregulated after drought, which likely drive the changes in cuticle composition and ultimately increased glyphosate resistance following drought. Ultimately, understanding how drought impacts glyphosate resistance is critical for maintaining optimal weed control in the changing climate. HighlightDrought stress induces changes to cuticle composition and gene expression that reduce glyphosate absorption, thereby increasing horseweeds ability to survive glyphosate application.

The fungal pathogen Pyricularia oryzae is notorious for causing blast disease in various important cereal crops, including wheat, rice, millet, and oat. Whole-genome-informed data on this pathogen are necessary to better understand the host adaptability of the fungus, including identifying key determinants of infection to enable more precise disease control. Here, we report highly contiguous genome sequences (using long-read PacBio technology) of two isolates from rice paddies in Portugal, M22.7 and T22.2, which exhibit distinctly aggressive symptoms in rice. Both mitochondrial and nuclear sequences were characterised in this study. The resulting nuclear genomes have assembly lengths of 46.4 Mb for M22.7 (198x coverage) and 46.3 Mb for T22.2 (163x coverage), with near-complete BUSCO completeness (98.8%) and a 0% contamination score (EukCC). Phenotypic analysis showed M22.7 to be more virulent than T22.2, which may be explained by the lower number of predicted effector genes and higher transposable element content in M22.7 relative to T22.2. This announcement represents the first genome resource for natural isolates of P. oryzae from Portugal in over 20 years, filling an important data gap from a major European rice-producing country that produces locally adapted rice varieties under specific agro-environmental conditions (near the Atlantic coast).

A highly contiguous and complete reference genome of Cercospora cf. flagellaris, the causal agent of foliar disease on many plant hosts including Cercospora leaf blight of soybean, was assembled using a combination of PacBio and Illumina sequencing reads. The genome assembly is 33.72 Mb in length and consists of 14 nuclear scaffolds and one mitochondrial contig. Four scaffolds have telomeric repeats on both ends and represent fully assembled chromosomes, while nine scaffolds represent partially assembled chromosomes with telomeric repeats on one end. The assembly has an N50 of 2.90 Mb and an L50 of 5 scaffolds. Genome annotation identified 11,268 genes, of which 947 and 360 were predicted to encode secreted proteins and effectors, respectively. Additionally, 512 genes were predicted to encode carbohydrate-active enzymes and 60 biosynthetic gene clusters were annotated. Taken together, this annotated genome assembly will be a valuable resource for genomics, host-pathogen interactions, and population biology research in this economically important pathosystem.

We present Nemo 2.4, an advanced forward-time individual-based simulation framework designed to model the complex eco-evolutionary dynamics and genetic basis of quantitative traits. This tool addresses current challenges in evolutionary quantitative genetics by providing unprecedented flexibility and computational efficiency. Nemo 2.4's modular architecture allows researchers to design custom life cycles by combining specialized Life Cycle Event (LCE) modules, from reproduction and dispersal to selection, crossing, and phenotype expression. The software supports diverse population models, including both Wright-Fisher (WF) and non-WF dynamics, spatially explicit models, and varying demography. Nemo 2.4 handles a wide range of genetic architectures, including both multi-allelic Quantitative Trait Loci (QTL) for general trait studies, and dense di-allelic Quantitative Trait Nucleotides (QTN) implemented with highly optimized bit-wise data structures. Crucially, it allows the simulation of QTNs on comprehensive genetic maps that incorporate other genetic elements, providing genomic-scale resolution. Key biological complexities are integrated natively: the model accommodates modular pleiotropy, dominance, and pairwise epistasis across multiple traits, facilitating the study of complex genotype-phenotype mappings. Furthermore, Nemo 2.4 models phenotypic plasticity through reaction norms and incorporates underlying liability thresholds, enabling the simulation of environmental influences on trait evolution with various forms of selection (e.g., Gaussian, linear, truncation). Due to its compiled design and memory-efficient data representations for large numbers of loci, Nemo provides a robust platform for running high-throughput simulations critical for testing theoretical predictions in polygenic adaptation and understanding evolutionary responses to changing environments.

Modern hybrid maize (Zea mays L.) breeding programs are based on the management of distinct complementary heterotic groups to maximize heterosis in high-performing hybrids. This practice lowers shared genetic segments and increases divergence between groups to limit inbreeding in hybrids. However, most breeding programs have not always enforced strict separation between heterotic groups in the past. Competitor commercial hybrids were notably a common elite germplasm source for inbred development, which would diminish divergence between groups. This study proposes a new haplotype-based approach to assess hybrids residual inbreeding based on parental similarity. The new haplotype method has a stronger significant negative effect on hybrids grain yield than raw SNP data. Evaluation of modern experimental hybrids uncovered related inbreds contributing to superior rates of residual inbreeding. Analysis of these inbreds revealed haplotype transfers between heterotic groups, originating notably from the use of a Stiff Stalk-Iodent commercial hybrid as breeding starts material in both Stiff Stalk and Non-Stiff Stalk breeding populations. The introduction of this intergroup parent generated heterotic-group-specific haplotype migration between crossing pools. These fragments caused significant genome-wide residual inbreeding in experimental hybrids across selection cycles. This study highlights the necessity for accurate evaluation of external sources of diversity to minimize haplotype transfers and admixture between crossing pools. We demonstrate the consequences of using commercial hybrids in inbred development, particularly regarding residual inbreeding, and their effects on hybrid performance. Insights from these results can assist breeders in optimizing the choice of parents for introducing genetic diversity in a reciprocal recurrent selection scheme. KEY MESSAGEHaplotype-based hybrids parental similarity better predicts grain yield than marker-based identity-by-state. Utilization of commercial hybrids as breeding start material resulted in higher hybrid residual inbreeding even after several selection cycles

Caenorhabditis nematodes are a powerful model clade for evolutionary genetics. Isogenic lines and panels of recombinant inbred lines (RILs) are among the most essential tools for genetic studies in these species. While most Caenorhabditis species are gonochoristic, large RIL panels have only been developed for self-fertilizing species. This gap biases our understanding and limits our ability to address questions related to the genetic architecture of traits in outbred populations, which have radically higher genetic diversity, heterozygosity, and effective recombination than selfers. Having previously identified Caenorhabditis becei as a tractable gonochoristic species due to its moderate inbreeding depression, we generated two panels of advanced-intercross RILs derived from three individual outbred C. becei worms collected from a single locality on Barro Colorado Island, Panama. One panel derives from a pair of worms sampled from a single rotting fig; the other derives from a cross between worms from two different figs. The panels share one founder in common, yielding two half-sib RIL panels. We sequenced and haplotyped the lines, identifying millions of variants and thousands of recombination breakpoints. Using simulations, we demonstrate the suitability of these lines for quantitative genetics studies and QTL mapping. In our single-fig panel, we observe abundant heritable variation in population growth rate, individual body size, and sexual dimorphism for body size. We detected four QTLs associated with population growth rate and show that estimated allelic effects are good predictors of selection that occurred during panel derivation.

Organisms have spent their life histories exposed to background levels of natural ionizing radiation. To document the role that radiation plays, the deprivation of these natural levels has been studied by incubating organisms in the shielded space of underground laboratories. We report here on two studies (Study I and Study II) using Aedes aegypti for the first time as a model organism incubated 655 meters underground at the Waste Isolation Pilot Plant (WIPP) outside of Carlsbad, New Mexico, U.S.A. Male mosquitos were incubated at the surface exposed to natural background radiation, and were compared to two underground treatments in which incubators were supplemented with radiation sources used to mimic background and these groups were compared to the underground, radiation-deprived treatment. In Study I, the mosquitos incubated underground in the absence of natural radiation had higher levels of mortality compared to those incubated at the surface and PCA plots of the two transcriptomes were clearly differentiated. Study II was conducted the following year and the experiment was narrowed to include only the surface control and underground, radiation-deprived treatment which allowed for four biological replicates. Again, there was a higher level of mortality in the mosquitos grown underground compared mosquitos grown at the surface. Transcriptomes were not as clearly differentiated by PCA analysis and fecundity data were similar between the two groups. Functional analysis of transcriptomic DEGs from two independent studies suggested there are stress responses in radiation deprived mosquitoes. The absence of a secondary stressor in Study II is discussed as an explanation for the transcriptome differences in the two experiments.

HighlightThe sweet fragrance of lettuce was attributed, for the first time, to the synthesis of 2-acetyl-1-pyrroline caused by a deficiency in the betaine aldehyde dehydrogenase gene. Fragrance is among the most valuable traits of high-quality crops and influences consumer preferences. Although 2-acetyl-1-pyrroline (2AP) is a key component of fragrant cultivars in several crops, its genetic mechanism in lettuce (Lactuca sativa L.) remains poorly understood. The betaine aldehyde dehydrogenase (BADH) gene has been identified as causative for 2AP-derived fragrance in rice and soybean cultivars. Hence, we conducted a linkage analysis using an F2 population derived from a cross between Kukichisya (fragrant) and Rennet (non-fragrant) for three candidate genes of BADH orthologs in the lettuce genome. Analysis linked LOC111877932 located in LG4 to the fragrance trait, and it was designated LsBADH1. Comparison among Kukichisya, Salinas, and candidate BADH of sunflower (Helianthus annuus L.) revealed three non-synonymous single-nucleotide polymorphisms (nsSNPs) in exons 1, 2, and 9, and suggested that nsSNP in exon 9 was strongly correlated with fragrance in Kukichisya. A premature stop codon introduced in exon 5 of LsBADH1 using Target-AID base-editing technology resulted in truncated BADH1 and higher 2AP levels. Our results indicated that LsBADH1 is responsible for the 2AP-derived fragrance. Our findings can be applied to select cultivars based on a novel concept for the cooking process, providing a transformative platform to breed fragrant lettuce as a high-value-added product.

Redox balance is central to aerobic metabolism, yet acute heat stress can destabilize this balance by increasing metabolic rates and shifting the balance of critical electron carriers such as NADH. In early Drosophila melanogaster embryos, maintaining redox balance is particularly critical as embryos undergo a developmental redox shift and rely on oxidative phosphorylation to power nuclear divisions. Here, we assayed six isofemale D. melanogaster lines from temperate (Vermont, USA; France; Japan) and tropical (St. Kitts; Ghana; India) climates to assess metabolic responses to heat in heat-sensitive versus heat-tolerant embryos. We used untargeted LC--MS to measure 33 metabolites and the major redox couples (NADH/NAD+, NADPH/NADP+, and GSH/GSSG) at 25{degrees}C and after a 32{degrees}C heat shock. In all embryos, heat shock induced shared shifts in metabolic profiles, with increases in nucleotide monophosphates (e.g., AMP, CMP, and GMP) and amino acids (e.g., alanine, glutamic acid, serine). In contrast, redox metabolites diverged by region: heat-sensitive temperate embryos shifted toward a more oxidized state (46.6% decrease in NADH/NAD+ ratio and 4-fold increase in oxidized glutathione), while heat-tolerant tropical embryos maintained glutathione balance and increased the NADH/NAD+ ratio by 52.9%, indicating a more reduced state. These patterns are consistent with higher NADH oxidation and greater oxidative stress (inferred from oxidized glutathione) in the temperate embryos, versus better maintenance of redox balance in tropical embryos. Together, our results suggest that maintaining redox balance is a key determinant of acute heat tolerance, and healthy development overall, during early embryogenesis.